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A general model for opinion formation and competition, like in ideological struggles is formulated. 
The underlying set is a closed one, like a country but in which the population size is variable in 
time. Several ideologies compete to increase their number of adepts. Such followers can be either 
converted from one ideology to another or become followers of an ideology though being previously 
ideologically-free. A reverse process is also allowed. We consider two kinds of conversion: unitary 
conversion, e.g. by means of mass communication tools, or binary conversion, e.g. by means of 
interactions between people. It is found that the steady state,when it exists, depends on the number 
of ideologies. Moreover when the number of ideologies increases some tension arises between them. 
This tension can change in the course of time. We propose to measure the ideology tensions through 
an appropriately defined scale index. 

PACS numbers: 



I. INTRODUCTION 
A. Nonlinear systems, opinion formation and population dynamics 

The case of opinion formation [I] - [5] in sociology is seen in recent physics works as analogous to state phase 
evolution in non-equilibrium systems. Applications to population dynamics [B], extinction of populations [7], animal 
and human migration |8], policy and politics p], languages |T0], religions [TTl - are found to be similar to 
epidemics [15 , forest fires [THj and other self-organizing systems much studied in statistical physics. No general 
pattern is however available and situations can be very varied. Whence it seems appropriate to continue considering 
the questions through the methods used in the theory of turbulence |17j - |20j . low-dimensional dynamical systems \2T\ 
- [23] or theory of nonlinear waves [21] - [26] but exordium to go back to the classical Verhulst ideas and Lotka-Volterra 
model by introducing some realistic conditions on the growth ratios and on the interaction coefficients between the 
populations [27] - |31j . Indeed there is then a connexion to the problem of extinction of populations religions [If] , 
languages |32j . and to the very modern question of internet governance in which the old stakeholders, i.e. the most 
powerful actors, and a variable set of new participants are somewhat abused or lacking cohesion in their reaction [33J. 

Interestingly this demands at some stage a consideration of the connexion between economical and social physics 
and social dynamics, including the analysis of time series [M] - [39] though with some caveat due to either size or 
debatable data value. Whence the need for a theoretical approach and some modelization in order to focus any data 
gathering toward useful input in further work. 

B. Organization of the paper 

In Sect. 2, we formulate a general/mathematical model in a finite system size, we emphasize, allowing for a changing 
population size. In a realistic way we consider people unaffected by the available ideologies as well as conversions to 
ideologies. Two mechanisms are discussed by which the followers of an ideology can increase: unitary conversion (a 
citizen is converted by means of forms of mass communication such as newspapers, radio or television channels) and 
binary conversion (a citizen is converted by interpersonal contacts with other citizens). 
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In Sect. 3 we study the case when only one ideology spreads among the population. The whole country population 
is found to evolve towards an equilibrium state. In this state some part of the people becomes followers of the ideology 
and the remaining ones are not followers of the ideology. The fraction of adepts depends on the intensity of the unitary 
and binary conversion as well as on the ability of the ideology to reduce the dissatisfaction among its followers. In 
Sect. 4, the case of two ideologies is examined. The introduction of a second ideology leads to tensions between 
the ideologies as the number of followers drops in comparison to the case when each of the ideologies is alone in the 
country. The tension can be quantified on a scale that can be considered to be an index of ideological tensions in the 
society. 

In Sect. 5, in the case of three ideologies, we briefly show that a chaotic behavior is seen to occur among many 
other solutions. This seems close to intuition. In Sect. 6, conclusions are outlined and connected for relevance sake 
to a usual observation that the coexistence of ideologies imphes the existence of tensions between them. 



II. MATHEMATICAL FORMULATION OF THE MODEL 



Let us consider a set ( country ) with a population of N agents. We are going to consider that the population is 
divided into n+1 factions: n factions each with a different specific ideology, such that the number of members in the 
corresponding populations are TVi, A^2, ■ • ■ , Nn, and a fraction A^o of people which are not followers of any ideology at 
a given moment of time. Then 

n 

N^No + Y,Ni (1) 
1=1 

We assume that the overall population evolves according to the generalized Vcrhulst law 

dN 



r{t,N,Ni,...,N„,p^)N X 
N 

^~ C{t,N,Ni,...,N„,Pf,) 



(2) 



where = (pi, . . . £^re parameters, describing the environment. 

The growth process is constantly disrupted by small extinction events, as in |3D], monitored through r{t,N,Ni, 
. . . ,Nn,p^). r is the overall population growth rate; r can be positive or negative. In this paper we shall consider 
the case r > 0, i.e. we shall study the spreading and competition between ideologies in a country with a growing 
total population. C'(t, N, Ni, . . . ,Nn,Pfj,) is the maximum possible population of the country (its so called carrying 
capacity). In every ideological population i we have to account for the following processes: deaths, dissatisfaction, 
unitary conversion, and binary conversion. 

1. First, we expect a decrease of the number of followers of an ideology through death or dissatisfaction with the 
ideology, i.e. through a term ViNi, where < 0. In general — riit, N, Ni, . . . , Nn,p^). 

2. Unitary conversion: such a conversion from one ideology to another is made without direct contact between the 
followers of different ideologies. The conversion happens through the information environment of the population. 
Elements of this environment for example are the newspapers, the radio stations, television channels, printed 
propaganda materials or mass events such as speeches during election campaigns. Excluded are the direct 
interpersonal contacts which lead to the binary conversion described below. In order to model the unitary 
conversion we assume that the number of people converted from ideology j to ideology i is proportional to 
the number Nj of the followers of the ideology j. An fij coefficient characterizes the intensity with which this 
conversion occurs. The corresponding modeling term is fijNj. We assume that fa = 0. In addition a term 
fioNf) describes the unitary conversion toward ideology i from the Nq people who were not followers of any 
ideology at the corresponding moment of time. In general 

/,j = /,,(t,iV,iVi,...,7V„,p^,C) 



f^o^ Mt,N,N,,...,Nn,p^,C) 

3. Binary conversion: bijkNjNk- In general this term describes the conversion to the i-th ideology because of direct 
interaction between members of the j-th and fc-th ideology. We assume that the intensity of the interpersonal 
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contacts is proportional to the numbers Nj and Nk of the followers of the two ideologies. The coefficient 
characterizing the intensity of the binary conversion is bijk- The larger is bijk, the more people are converted 
to the i-th ideology. In general the binary conversion coefficients can be bijk = bijk{t, N,Ni, . . . , NniP^, C). Of 
course bm = 0: there is no self-conversion. 

The equation for the evolution of the followers of the ideology i becomes 

'^=ri{t,N,Nu...,Nn,p^,C)Ni + 
Mt,N,Ni,...,Nn,p^„C)No + 

n 

J2fi,{t,N,N,,...,Nn,p^,C)N,+ 
bio{t, N,Nu..., Nn,Pi„ C)NiNo + 

n n 

^ijk{t, N,Ni,..., 7V„,p^, C)NjNk (3) 

j=i fe=i 

In general we can have a co-evolution of the environment and the populations, i.e. = Pij,{N, Ni,. . . , N^, C, t), but 
this will not be discussed here. 

Indeed in this paper we shall discuss the simplest version of the model namely the case in which all the coefficients 
are time and p^ independent. Then the model system becomes geared by 



N = No + J2^i (4) 



i=l 



= nNi + fioNo + fijNj + hoNiNo + ^ijkNjNk (6) 

j=i j=i fe=i 

Two remarks are in order. 

1. Notice that arbitrary values are not allowed for the coefficients of the model. They must have values such that 

iV, iVo, A^i, . . . , Nn be nonnegative at each t. 

2. Let us consider the i-th population and the binary conversion characterized by the coefficients bijk and bikj 
where j and k are different from i. One could at first think that bijkNjNk and bikjN^Nj describe one and the 
same process in which the interaction between followers of the j'-th and fc-th ideology leads to a conversion of 
these to the i-th ideology. In general however one should not identify the two terms. In so doing in the general 
model we retain one additional degree of freedom, i.e., that which allows to distinguish between the ideology 
that is of the initiator of the interaction and the ideology of someone who is apparently simply a participant in 
the interaction. 

Below we shall consider the dynamics of populations of followers of the ideologies for the cases of presence of 1,2 or 
3 ideologies in the country. 



III. THE CASE OF ONE IDEOLOGY 



In the case of spreading of one ideology the population of the country is divided into two groups: Ni followers of 

the single ideology and people who are not followers of this ideology at the corresponding time. Let us first discuss 
the case when only the unitary conversion scheme exists, as possibly moving the ideology-free population toward the 
single ideology, i.e. /lo is finite. Let the initial conditions be N{t = 0) = N{Q) and Ni{t = 0) = A''i(0). The solution 
of the model system is 

^ ' 7V(0) + (C - iV(0))e-'-* ^ ' 
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Ni{t) = e-(/i«-''i)*<{ Ni{0) + ^ X 



'J C~N{0) /lo - ri 
7V(0)e'-* ' ' r 



(8) 



where $ is the special function 



$(z,a,w) 



Y— 

^ (v + 

n— 



Z |< 1 



(9) 



(10) 



The condition | z |< 1 is equivalent to N* > C/2. This means that the solution (7) describes the development in the 
case of a densely populated territory. For N* < C/2, $ can tend to oo and we can only obtain a numerical solution 
of the model system of equations. 

The obtained solution describes an evolution in which the total population N reaches asymptotically the carrying 
capacity C of the environment. The number of adepts of the ideology reaches an equilibrium value which corresponds 
to the fixed point of the model equation for This fixed point is 



Ni = 



C/io 

fio - ri 



The number of people which are not followers of the ideology asymptotically tends to A'o = C — iVi . As a numerical 
example let C — 1, fio = 0.03 and ri — —0.02, then Ni = 0.6 which means that the evolution of the system leads to 
an asymptotic state in which 60 % of the population are followers of the ideology and 40 % are not. 

Now let not only unitary but also binary conversion processes be possible. The evolution in this case cannot be 
investigated analytically. However the asymptotic behavior for A'^i can be obtained when the total population N has 
reached the carrying capacity C of the environment. For this asymptotic state the evolution of Ni reads 



dNi 



riiVi + /io(C - iVi) + hoNiiC - Ni) 



There exist two fixed points but only one of them satisfies the requirement iVi > 0. This fixed point is 

in - fio + hoC) + V(?-i-/io + 6ioC)2 + 46io/ioC 



26i 



(11) 



(12) 



The equation (11 1 has an analytical solution. The solution depends on whether Ni > Ni or A^i < A^i. The two cases 



can be reahzed respectively when iVi(O) > A^i and Ni{0) < Ni. If A^i(O) > Ni then 



(13) 



where 



Xi^n- fio + bwC + V(^i - /lo + bioCy + 46io/ioC + 

g-(i+r)(ri-/io+6ioC) ^ 

[ri - fio + bwC - - /lo + bioCy + 46io/ioC^ 



Yi = 2b 
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(a) (b) 




250 500 750 1000 250 500 750 lOOO 

t t 

FIG. 1: (a): Illustration of the inertial growth and its dependence on the parameter n. The other parameters are: C = 1, 
r — 0.2, /lo — hio = 0.001. (b): Illustration of the inertial growth and its dependence on the parameter r, with C = 1, 
n = -0.02, /lo = 6io = 0.01 



and where the characteristic time r is here given by 



In 



1 

n - /lo + hoC 
2biojVi(0) + Z 
26ioiVi(0) + ri-/io + 6ioC-Z 



(14) 



where Z = ^y {n - fw + hwCy + 46io/ioC. 
For the case < Ni{Q) < Ni 



iVi = # (15) 
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where 



X2 — Xi 



Y2 = 2ho [1 + e-('+^>in - fio + hoC) 

Let us now discuss the time behavior of the number of followers of the ideology. There are three possibilities 

1. > for all t, i.e. the number of followers increases monotonically. For the particular case of only a unitary 
conversion process the condition reads < ^ 

2. < for all t, i.e. the number of followers decreases monotonically. For the case of only unitary conversions, 
the condition reads — -P^ > ^ 

710 JVi 

3. The most interesting case is when can change sign with increasing t. The following effect can be observed: 
the number of followers of the ideology can increase despite the fact that ri < 0. The reason for this effect is 
the increasing number A^O) occurring because of the fast enough growth of the population of the country. When 
A^o is small the term containing ri dominates and A^i decreases. But in the c;oursc of time A^o increases. Then 
the conversion begins to dominate over dissatisfaction and the number of the followers of the ideology begins to 
increase. We shall call such a kind of growth of the followers of the ideology an inertial growth. 

Fig. 1 illustrates an inertial growth. In Fig. la, it can be observed that the process of initial shrinking and then of 

inertial growth can exist for a large range of coefficient values. Inertial growth can exist even if the ideology is weak 
with respect to its keeping of followers, i.e. when n has large negative values. Note that the figure illustrates the 
case when the carrying capacity of the environment is a constant. If the carrying capacity would change, one could 
observe sequences of phases of inertial growth and shrinking. Fig. lb shows that a small growth rate r, i.e. of the 
total population, leads to a slowing down of the inertial growth process. 
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IV. CASE OF TWO IDEOLOGIES: THE IDEOLOGICAL TENSION 

In this section we will discuss the competition for adepts that the presence of a second ideology introduces; this is 

leading to a measurable conflictual tension, as will be shown. 

Let us consider the model system for the case of two ideologies with populations of followers A^i and N2; we assume 
that all parameters are kept constant. We have 

at 

(6112 + 6l2l)iVliV2 + 6l22A^| (17) 

^ = r2iV2 + /2oA^o + /2iA^i + h2oNoN2 + 

b2nN^ + {b2i2 + b22i)NiN2 (18) 



N = N0 + N1 + N2 (19) 

Let us quantify the tension between ideologies by means of an asymptotic analysis. We discuss the case of only unitary 
conversion of members of the population that are not followers of any of the ideologies. In order to emphasize the 
unitary conversion effects let us assume that the binary conversion as well as the unitary conversion from one ideology 
to the the other one are negligible. In this case /12 = /21 = and 6112 = 6122 = &121 = ^211 = &212 = ^221 = 0. In 
addition let us consider the asymptotic case in which the total population has reached the carrying capacity : N = C. 
The equilibrium state is characterized by the fixed points : 

m= , f;^" — , ^2^ , f/" — (20) 

If the ideologies were without competition their size would be; see the previous section 

N, = ^^; N2 = -^ (21) 
/lo - r\ /20 - r-2 

It can thus be observed that the popularity of an ideology shrinks when competing ideology or ideologies spread 
around the country. Let us evaluate this shrinking. We have 



TVi 1 iV; 



2 



1 J- ''1/20 ' AT 1 I ''2/10 

^ + r2(/io-ri) ^2 1+ r-i(/20-r2) 



(22) 



As a numerical example let ri = r2 = —0.01 and /lo = /20 = 0.02. Then Ni/Ni = 0.6, i.e. the number of followers 
of the ideology 1 descreases by 40%. Of course this causes some tension between the ideologies. A measure of this 
tension can be through the index 

Ti,k = l-^, (23) 

(k) 

where is the population of the followers of the i-th ideology when the fc-th ideology is presented in the country 

too. If the ideology is alone then n[^^ = Ni and the tension index is Ti-i — 0. If Ni decreases because of the 
competition with the second ideology, then the tension between the ideologies characterised by the tension index T^k 
increases. The above definition for the tension holds even if A''! follows some time dependent trajectory. 

The tension index can be generalized for the case of an arbitrary number of ideologies in the country (next section 
for example). Let m ideologies be presented in the country. The tension on i-th ideology in presence of two other 
ideologies, k and I is 

T,.,k,i{t) = l-^^ (24) 
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FIG. 2; Evolution of the total population A'^, number of followers of the ideologies, A'^i and and the number of A'^o of people 
that are not followers of any ideology, (a): ri = -0.03. Initial conditions: Afi(O) = 0.03, iV2(0) = 0.012, iV(0) = 0.2; (b): 
ri = -0.02. Initial conditions: iVi(O) = 0.3, iV2(0) = 0.18, iV(0) = 0.5. In both figures : = -0.005, ra = 0.01, C = 1, 
/lo = 0.001, /20 = 0.003, /ii = /22 = 0, /12 = /21 = 0.001, 610 = &20 = 0.001, bm = 6222 = 0; 6112 = &121 = &122 = &211 = 
&212 = fe221 = 0.001. 



where ^^^'^^ is the population of followers of the i-th ideology when the ideologies fc and I operate in the country too. 
The tension on i-th ideology in presence of three other ideologies, j, k, I is 



{3,k.l) 



it) 



N,, 



(25) 



where iV^^^*^''-* is the number of the followers of the i-th ideology in presence of ideologies j, k, I. In such a way we 
can define a series of indices for the quantification of the tensions among ideologies competing for followers ( in the 
same country). 

Fig. 2 shows typical results from the numerical investigation of ( 16 ) -( 19 ). Fig. 2a shows the purely inertial growth 



of a population of followers of ideology 2 and its decline, followed by the inertial growth of the population of followers 
of ideology 1. In Fig. 2b one can observe an initial decline followed by an inertial growth of the number of followers 
of both ideologies. These results can be usefully compared to a model of telecommunication competition |45] where 
it is concluded that (we quote) "schemes targeting local cliques within the network are more successful at gaining a 
larger share of the population than those that target users randomly at a global scale (e.g., television commercials, 
print ads, etc.). This suggests that success in the competition is dependent not only on the number of individuals in 
the population but also on how they are connected in the network. The network in our above investigation is such 
that all agents are fully connected with each other; we consider a fully connected graph. This is equivalent to a mean 
field approximation study. Notice that the links are weighted through the fij and bijk coefficients. 



CASE OF THREE IDEOLOGIES 



The case of three ideologies will not be treated in full here. It can be easily understood that the above procedure 
can be extended to the case of multiple ideologies. We only restrict our presentation to the main points and illustrate 
the newness. In the case of three ideologies the model system Q - (|6]) has an analytical solution in the following 
conditions. First, let the binary conversion be negligible: bijk = 0, bio = 0. Then the model system reduces to 



(26) 



(27) 



dt 



= {n - no)N, + f.oN + 



(28) 
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FIG. 3: An example of chaotic behavior of the tension indices for the case of three competing ideologies. The model system 
is as in Eq.(32l. For the attractor shown the parameters are as follows: kh — ki,ki2 = ki3 ~ k,2, K2i ~ — ki, K22 = —k,2, 



K.23 ~ 1^2, K31 = K3, K32 = K2, K33 = K2. C = 10, Kl = 0.5, K2 = 0.1, = 1.63, 610 = 0.001, 620 = 0.002, 630 = 0.001, 
n = -0.01,r2 = -0.03,r3 = -0.01. 



Let the initial conditions be N{t = 0) = Af(0), iVj(t = 0) = iVj(O) and N{id) > C/2. Then the solution is 

CN{0) 



Nit) 



N{0) + {C~N{0))e' 



(29) 



c/io 



et(/.o-n)$ _ 



C-NjO) f,o - r, 

Af(0) ' ' r 

g-iV(O) ho - r-. 

iV(0)e'-* ' ' r 



(30) 



7Vo(i) = 7V(t)-^iV,(t) 



(31) 



where <i> is the special function defined in section 3. Above i = 1,2, ... ,n. When n ~ I one returns to the case of 
section 3. For n — 2 one recovers the solution of the case of section 4. 

Chaotic change of the ideological tensions becomes possible when the number of the ideologies becomes equal to 3 
or larger. This can be demonstrated easily as follows. Let the total population having reached the carrying capacity 
of the environment N = C. Let /12 = /ra = /lo — and bijk — b2jk — b^jk = for j, k — 1,2, 3. Let also 
/21 = /23 = /20 = ™d /31 = /32 = /30 = 0. Let us rewrite the parameters in the following way 

ri + /ii Kii + Kr2 + K13, ^111 - &10 -Kii 
bii2 ~ ^10 — —K12, bii3 — biQ — — K13 

r2 + /22 = K21 + K22 + ^23,^222 - ^^20 = ~ ^22 



&22r — ^20 — —^21, &223 " ^30 — —^23 



7'3 + /33 = '«31 + K32 + K33, 6331 - &30 = -K3I 



&332 — ^30 — ~'*32j ^333 ~ ^30 — ~>^33 
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The model system of equations becomes 



^ = 7V,^K,,(l-iV,) + 6,oC^^»; i,j = 1,2,3. (32) 



The existence of chaos in a particular case of ( 32 1 was discussed in [46 . Fig. 3 is an illustration of a case of chaotic 
change of the ideological tensions for the (32 1 corresponding model. is always significantly different from 0. This 

means that because of the competition with the ideologies 2 and 3 the ideology 1 remains at a significant " distance" 
from its most favorable state, - this one which would exist if which there was no competitor. In addition the ideology 
1 copes relatively good with the situation as its tension index remains relatively distant from 1. The other two 
ideologies experience large tensions and from time to time are close to extinction. The ideology 2 evolves better than 
the ideology 3 which experiences large oscillations of the number of followers as consequence of the competition with 
the other ideologies. This illustration for a given set of parameters indicates the interest of the approach, since a 
few measures would allow to calibrate the parameters, in specific situations, whence would lead to considerations 
pertaining to forecasting science. 



VI. CONCLUDING REMARKS 



In this paper a general model for ideological competition was formulated. The model applies to cases in which 
countries have a variable total population which evolves according to the generalized Verhulst model. The discussion 
in the paper is concentrated on the cases of constant coefficients and when the total population of the country 
increases. An original ingredient concerns also the number of followers of an ideology which can increase without 
interpersonal contacts, but solely on the basis of so called unitary conversion, e.g. as a result of different forms of mass 
communication. The number of followers can also increase by means of binary conversion as a result of interpersonal 
contacts. It is emphasized that the conversion can be outside the competing ideologies of interacting agents. 

The dynamics of the populations of followers of the ideologies is discussed for the case of one ideology and for the 
case of two and three competing ideologies, in Sect. 3-5, respectively. For the case of one ideology, the simple version 
of the general model describes the evolution to an equilibrium state in which the population consists of some amount 
of followers of the ideology and persons indifferent to the ideology. The introduction of a second ideology leads to 
some tension between the ideologies as the numbers of followers drop in comparison to the case when each of the 
ideologies is alone in the country. The ideological tensions can be quantified by a set of indices. A nonzero index is 
a characteristic feature of the competition. Each ideology most of the time, if not always, tries to set its indices of 
tension to 0, i.e. it tries to reach its maximum number of followers (which is the case when the ideology is alone in 
the country) . This can be done by decreasing the number of followers of the other ideology on the territory. 

We have indicated that chaos can exist (can be found) when the number of available ideologies increases above 
2. The number of parameters is considerable, as in many realistic population evolution studies. However the set 
of parameters appears to be realistic enough to be calibrated in specific situations. This would lead to forecasting 
considerations. This and the other obtained results hint to good perspectives for applications of the methods of 
statistical physics, theory of networks, sociophysics, etc [HJ to the problems of ideological competition. This 

will be a subject of future research. 



Acknowledgments 



We would like to thank the ESF Action COST MP0801 (Physics of Competition and Conflict) for support of 
our research. Special thanks are devoted to Noemi Olivera, Andrew Roach, Jiirgen Mimkes and John Hayward for 
stimulating discussions on the dynamics of ideological competition. 



[1] S. Galam, F. Jacobs. Physica A 381, (2007) 366-376. 

[2] B. Pabjan, A. Pekalski. Physica A 387, (2008) 6183-6189. 

[3] S. Galam. Int. J. Mod. Phys. C 19, (2008) 409-440. 

[4] K. Sznajd-Weron, J. Sznajd. Int. J. Mod. Phys. C, 11, (2000) 1157 - 1165. 

[5] J. A. Holyst, K. Kacperski, F. Schweitzer. Physica A 285, (2000) 199-210. 

[6] M. P. Hassel, R. M. May. Nature 353, (1991) 255-258. 



10 



[7] A. Pekalski, M. Ausloos. Physica A 387, (2008) 2526-2534. 

[8] N. K. Vitanov, I. P. Jordanov , Z. I. Dimitrova. Communications in Nonlinear Science and Numerical Simulations 14, 

(2009) 2379-2388. 
[9] S. Moss. Proc. Natl. Acad. Sci. USA 99, (2002) 7267-7274 
[10] L. DallAsta ,A. Baronchelli, A. Barrat, V. Loreto . Non-equilibrium dynamics of language games on complex networks 

http : // arxiv .org /pdf /physics /0607054: 
[11] M. Ausloos, F. Petroni. Eur. Phys. Lett. 77, (2007) 38002 
[12] S. Picoli, R. S. Mendes. Phys. Rev. E 77, (2008) 036105. 
[13] J. Hayward. J. Math. Sociology 23, (1999) 255-292. 
[14] J. Hayward. J. Math. Sociology 29, (2005) 177-207. 

[15] R. M. Anderson, R. M. May. Infectious Diseases of Humans: Dynamics and Control. Oxford University Press, Oxford 
(1992). 

[16] J. A. M. S. Duarte. in Annual Review of Computational Physics, D. Stauffer, Ed. V, (1997) 1-23. 

[17] P. Berge, Y. Pomeau, C. Vidal. Order within Chaos: Towards a Deterministic Approach to Turbulence. Wiley, New York 
(1984). 

[18] N. P. Hoffmann, N. K. Vitanov. Phys. Lett. A 255, (1999) 277-286. 
[19] N. K. Vitanov. Phys. Rev. E 61, (2000) 956-959. 

[20] R. Z. Sagdeev, G. M. Zaslavsky. Nonlinear Physics from Pendulum to Turbulence and Chaos. CRC Press, Boca Raton, FL 
(1988). 

[21] M. A. Novak, R. M. May. Nature 359, (1992) 826-829. 

[22] S. Panchev, T. Spassova, N. K. Vitanov. Chaos Solitons & Fractals 33, (2007) 1658-1671. 
[23] E. Ott. Chaos in dynamical systems. Cambridge University Press, Cambridge (2002). 

[24] R. K. Dodd, J. C. Eilbeck, J. D. Gibbon, H. C. Morris. Solitons and nonlinear wave equations. Academic Press, London 
(1982). 

[25] N. K. Vitanov. Proc. Roy Soc. London A 454, (1998) 2409-2423. 

[26] M. Remoissenet. Waves called Solitons. Springer, Berlin (2003). 

[27] Z. L Dimitrova, N. K. Vitanov. Phys. Lett. A 272, (2000) 368-380. 

[28] Z. L Dimitrova, N. K. Vitanov. Physica A 300, (2001) 91-115. 

[29] Z. L Dimitrova, N. K. Vitanov. J. Phys. A: Math. Gen. 34, (2001) 7459-7473. 

[30] Z. L Dimitrova, N. K. Vitanov. Theoretical Population Biology 66, (2004) 1-12. 

[31] N. K. Vitanov, Z. I. Dimitrova, H. Kantz. Phys. Lett. A 349, (2006) 350-355. 

[32] D. M. Abrams, S. H. Strogatz. Nature 424, (2003) 900. 

[33] L. N. Olivera, A. N. Proto, M. Ausloos. Modeling the information society as a complex system, in press (2009). 
[34] Y. Aruka, J. Mimkes. Evolutionary and Industrial Economics Review 2, (2006) 145-160. 

[35] B. K. Chakrabartri, A. Chakraborti, A. Chaterjee. Econophysics and sociophysics: Trends and perspectives. Wiley- VCH 

Verlag, Weinheim, Germany (2006). 
[36] N. K. Vitanov, K. Sakai, L P. Jordanov, S. Managi, K. Demura. Physica A 382, (2007) 330-335. 
[37] G. Rotundo, M. Ausloos. Physica A 373, (2007) 569-585. 

[38] N. K. Vitanov, K. Sakai, Z. I. Dimitrova. Chaos Solitons & Fractals 37, (2008) 187-202. 
[39] A. Fabretti, M. Ausloos. Int. J. Mod. Phys. C 16, (2005) 671-706. 
[40] C. Wilke, T. Martinetz. Phys. Rev. E 56, (1997) 7128-7131. 

[41] J. Miskiewicz, M. Ausloos. (2006) Influence of information flow in the formation of economic cycles, in The Logistic Map 

and the Route to Chaos, M. Ausloos and M. Dirickx, Eds. (Springer, Berlin) (2006), 223-238 
[42] J. Miskiewicz, M. Ausloos. Physica A 382, (2007) 179-186 
[43] A. Lipowski, K. Gontarek, M. Ausloos. Physica A 388, (2009) 1849-1856. 

[44] J. Mimkes. A Thermodynamic Formulation of Social Science in Econophysics and Sociophysics, B. K. Chakrabarti, A. 

Chakraborti, and A. Chatterjee, Eds. (Wiley- VCH, Berhn) (2006). pp. 279-310. 
[45] E. F. Legara, A. Longjas, R. Batac. Int. J. Mod. Phys. C 20, (2009) 1-7. 
[46] A. Arneodo, P. Collet, C. Tresser. Phys. Lett. A 79, (1980) 259-263. 



